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Abstract 


Tbe  electrochemical  chargiiis  of  Pd  rods  by  deuterium  involves  a  complex  coupling  of  electrochemical,  interfadal  and  transport 
processes.  In  order  to  predict  the  overpotential,  surface  coverage  and  bulk  loading  of  the  electrode  during  charging,  a  model  has 
been  developed  that  incorporates  the  essential  features  ct  these  processes  and  involves  variables  such  as  the  electrochemical  rate 
constants,  the  bulk  difftisioa  coefficient  and  the  charging  current.  Features  of  the  computed  time  dependence  of  the  bulk  loading 
are  then  compared  with  pubhahed  experimental  charging  curves.  New  microscopic  obseivations  and  X-ray  diffraction  data  provide 
further  evidence  for  the  details  of  the  chargiiig  process. 


1.  lotroduction 

In  an  eailier  paper  [1]  we  presented  a  model  de¬ 
scribing  the  chaining  of  Pd  rods  by  eiectrochemicaliy 
generated  deuterium.  This  model,  which  is  based  upon 
the  coupling  of  interfadal  processes  with  the  diffu- 
skmal  transport  of  interstitial  offers  certain  predictive 
capabilities  with  respect  to  the  charging  process.  In 
particular,  we  examined  the  effect  of  the  magnitude  of 
the  rate  constants  of  the  electrochemical  steps  and  the 
charging  current  on  the  surface  coverage,  the  electrode 
potential  and  the  electrode  charging.  Subsequently, 
Riley  et  al.  [2]  reported  tm  the  experimentally  o^rved 
absorption  deuterium  into  the  Pd  lattice  of  cathodi- 
cally  polarized  electrodes  and  presented  results  with 
which  the  predictions  of  our  model  are  qualitatively  in 
good  agreement.  We  take  this  agreement  as  evidei^e 
for  the  general  validity  of  the  model  and  its  assump¬ 
tions  and  use  the  results  as  a  basis  for  a  refin^ 
interpretation  of  the  diarging  process. 

While  our  previous  discussion  of  the  model  [1]  was 
essentially  limited  to  the  oounination  of  the  behawior 

a  set  of  coupled  differential  equations,  this  paper 
presents  the  compariscm  of  the  m^el  and  its  pre^- 
tions  with  experiment  using  the  data  of  ref.  2.  We 
explore  the  relationship  between  the  various  operating 


fluxes  and  the  initial  charging  rate  as  well  as  the 
asymptotic  level  of  the  electrode  charging.  We  include 
new  results  based  upon  an  examination  of  microscopic 
observations  and  X-ray  difh'action  data  of  a  charged 
Pd  electrode.  Finally,  we  conclude  with  a  discussion  of 
the  strengths  and  weaknesses  of  the  model  in  the  light 
of  the  experimental  evidence. 

2.  Interpretation  of  experimental  data  of  Riley  et  al. 

The  data  published  by  Riley  et  al.  [2]  show  four 
characteristic  features  of  the  D/Pd  system  during 
charging:  (i)  the  existence  of  an  apparent  threshold 
value  for  the  cathodic  current  density  (CD)  above 
which  the  initial  charging  rate  does  not  increase,  imply¬ 
ing  the  existence  of  a  CD-dependent  charging  mod^s); 
(ii)  an  unexpected  dependence  of  the  asymptotic  elec¬ 
trode  loading  on  the  charging  CD,  i.e.  the  appearance 
of  a  maximum  (ref.  2,  Fig.  S);  (iii)  a  slower  rate  for 
“unloading”  than  for  “loading”,  which  implies  that  the 
rate  of  transport  between  the  metal  electr^e  and  bulk 
electrolyte,  i.e.  across  the  interphase,  depends  on  the 
direction — an  unlikely  situation  for  diffiision  control; 
(iv)  an  initial  charging  rate  that  depends,  for  the  most 
part,  inversely  on  the  diameter  of  the  Pd  electrode. 
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In  their  interpretation,  by  us&uming  diffusion  con¬ 
trol,  Riley  et  ai.  have  considered  only  events  occurring 
in  the  electrode  interior  and  have  disregarded  the 
influence  of  a  non-autonomous  interphase  on  the 
charging  process.  In  what  follows,  we  di^uss  their  data 
in  the  context  of  our  model  and  offer  a  somewhat 
different  interpretation  from  that  provided  in  ref.  2, 
showing  that  (i)  and  (ii)  are  natural  consequences  of 
the  model  while  (iii)  and  possibly  (iv)  require  the 
presence  of  more  complex  mechanisms. 

In  the  interpretation  of  the  charging  data  of  ref.  2,  it 
is  convenient  to  introduce  the  concept  of  the  charging 
efficiency  for  a  cylindrical  electrode  of  radius  r,  de¬ 
fined  here  as  the  ratio  of  the  absorbed  deuterium  to 
that  which  is  generated  by  electrolysis  at  a  current 
density  j  within  a  time  interval  At  and  given  by 


Z„F  r  AQ 
2  j  At 


(1) 


where  Z„  is  the  maximum  number  of  available  sites 
per  unit  volume,  F  is  the  Faraday  constant  and  AQ  is 
the  corresponding  change  in  the  loading  ratio  Q  of 
deuterium  to  palladium,  conveniently  expressed  in 
terms  of  the  D/Pd  atomic  ratio  since  is  taken  here 
to  be  the  numter  density  of  Pd.  The  first  factor  on  the 
right-hand  side  in  eqn.  (1)  is  a  constant  (equal  to 
5.45  X  10^  C  cm"’  for  a  Pd  electrode)  independent  of 
the  manner  in  which  the  experiment  is  carried  out;  the 
second  factor  is  controlled  by  the  investigator.  The 
third  factor  is  the  measured  charging  rate  and  reflects 
the  time  dependence  of  the  charging  efficiency  e  for 
for  an  eiqTeriment.  The  initial  charging  efficiency  char¬ 
acterizes  the  interphase  processes  independently  of  the 
electrode  geometry,  so  that  the  initial  charging  rate 
varies  with  r  as  1/r.  This  relationship  is  independent 
of  the  details  of  the  model. 


Fig-  1.  Coooept  of  the  Pd/D20  interphaie  (14,7]:  •,  adsofptkm 
layer,  t,  cbarie  tniufer  layer,  e,  surfue  divkhna  the  interphase  and 


the  hoik  electrode;  I,  lattice;  J,  total  current;  Rli,  Heyiovslqr-iioriuti 


path;  T,  Tafel  path;  V,  Volmer  path:  ^  and  >_«,  adioibed-adsoihed 


eachanae;  A.  non-autononaom  interphase  region. 


table  1.  Input  data  for  model  calculation  (cf  ref.  9) 


Parameter 

Value 

10,01 

s  Sx  lO"’  moJern"* 

lODl 

I.Dx  10  '*  mol  cm"  ’ 

lD,(s)l 

8..^x  10^’  mol  cm"’ 

1000  cm’  mol" '  s“  ' 

*2 

100 cm’  mol"  ‘  s"  ' 

KXIOcm’  mol"'  s"  ‘ 

*« 

I  .Ox  lO’  cm’  mol "  ‘  s " ' 

z- 

0.1 1.1  mol  cm"’ 

F. 

1.4 X  10  ■ "  mol  cm"  ’ 

C 

4.0X  10"’  Fcm"^ 

a 

0..S 

0.3 

*0 

l.Ox  10"’ 

5.0X10"’ 

N 

10 

The  processes  associated  with  the  electrolytic  charg¬ 
ing  of  the  Pd  electrode  are  shown  in  Fig.  1.  As  indi¬ 
cated,  there  are  six  fluxes  in  operation,  resulting  in  a 
large  number  of  governing  parameters.  Because  of  the 
large  number  of  parameters  involved,  the  predicted 
results  are  not  necessarily  computed  with  the  same 
values  of  CD  or  electrode  radius  as  those  for  which  the 
experiments  were  done.  Rather,  a  limited  number  of 
values  are  chosen  to  illustrate  possible  model  behavior, 
with  the  idea  in  mind  that  a  more  eleborate  procedure 
for  choosing  the  model  parameters  could  produce  a 
better  experimental  fit  if  in  fact  the  model  does  exhibit 
the  proper  behavior. 

llie  model  parameters  that  we  have  used  here  are 
listed  in  Table  1.  Several  of  the  parameters  differ  from 
those  used  earlier  in  ref.  1  where  we  provided  a  quali¬ 
tative  description  of  Pd/D  charging.  In  this  paper,  in 
order  to  obtain  agreement  with  the  experimental  data 
reported  in  ref.  2,  we  have  chosen  model  parameters 
accordingly.  Of  particular  significance  are  the  changes 
in  r„,  00  (o-  surface  density  of  sites,  r„,  was 

taken  to  be  0.005  of  the  maximum  number  of  surface 
sites  consistent  with  our  choice  of  Z„.  The  initial 
values  of  surface  coverage  and  bulk  loading,  Bq  and  (q, 
were  chosen  to  be  consistent  with  the  observation  of 
Riley  et  al.  [2]  that  the  equUibrated  loading  is  essen¬ 
tially  zero. 

Zl.  ThreshtM  CD;  charging  efftciency;  asymptotic  charg¬ 
ing  level 

Using  the  parameters  in  Table  1  and  the  diffusion 
coefficient  for  /3-PdD,  the  results  of  model  calculations 
for  a  cylindrical  electrode  of  radius  0.05  on  charged  at 
current  densities  ranging  from  30  to  240  mA  cm"’  are 
shown  in  Figs.  2  and  3.  It  is  seen  that  the  calculated 
charging  curves  reproduce  well  the  experimental  data 
reported  in  ref.  2,  Fig.  5.  In  particular.  Fig.  2  is  a  family 
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Fig.  2.  Model  calculation  of  loading  curves  as  a  function  of  current 
density:  electrode  radius  O.QS  cm;  current  densities  are  indicated  in 
mA  cm”*;  diffusion  coefficient  1.6x10”*  cm*  s”',  other  parame¬ 
ters  are  listed  in  Table  1. 


of  charging  curves  as  a  function  of  CD  and  illustrates 
that  the  model  predicts  a  saturation  of  the  initial 
charging  rate.  Figure  3,  curve  A  plots  the  initial  charg¬ 
ing  efficiency  calculated  by  eqn.  (1)  and  Fig.  3,  curve  B 
the  asymptotic  loading  level  obtained  from  the  curves 
of  Fig.  2.  The  model  predicts  a  brief  period  of  rapidly 
varying  charging  rate,  depending  on  the  electrode  ca¬ 
pacitance,  which  would  not  be  observable  experimen¬ 
tally,  therefore  the  initial  charging  rates  used  in  eqn. 
(1)  were  obtained  graphically  from  expanded  versions 
of  the  data  shown  in  Fig.  2  and  consequently  contain  d 
measurement  error.  In  terms  of  the  model,  what 'has 
been  referred  to  in  ref.  2  as  a  threshold  effect  in  the 


9 

I 


Fig.  3.  Model  calculations  as  a  function  of  current  density:  a.  initial 
charging  efficiency:  h.  asymptotic  loading  level:  eleetnxie  radius  0.05 
cm;  diffusion  coefficient  l.6x  10”*  cm*  s  ':  other  parameters  are 
listed  in  Table  I . 


charging  rate  with  increasing  CD  may  in  reality  be  a 
saturation  of  the  charging  rate  caused  by  a  decrease  in 
the  effectiveness  of  the  Volmer  path  with  respect  to 
the  Heyrovsky-Horiuti  path  for  charge  transfer,  coinci¬ 
dentally  resulting  in  a  charging  efficiency  that  appears 
to  be  inversely  proportional  to  the  current  density. 
Further,  the  maximum  in  the  measured  loading  curve 
asymptotic  levels  as  a  function  of  CD,  i.e.  characteristic 
feature  (ii)  (see  Fig.  3  and  eqn.  (11)  of  ref.  1),  probably 
results  from  the  same  shift  in  the  charge  transfer 
rather  than  from  changes  in  the  Tafel  path  or  the 
absorption  of  deuterium,  since  only  the  Volmer  and 
Heyrovsky-Horiuti  paths  (eqns.  (1)  and  (11)  of  ref.  1 
respectively)  exhibit  a  CD  (overpotential)  dependence. 

It  is  of  interest  to  note  that  earlier  calculations  that 
did  not  include  desorption  (the  Heyrovsky-Horiuti 
path)  showed  an  instability  in  the  computed  overpoten¬ 
tial,  suggesting  that  the  shift  in  balance  between  ad¬ 
sorption  and  desorption  plays  an  essential  role  in  the 
electrode  charging.  However,  a  direct  computation  of 
the  charging  currents  for  the  two  paths,  which  we  have 
not  done,  would  be  desirable  to  clarify  this  behavior. 

2.2.  Electrode  loading  vs.  unloading 

Evidence  of  an  active  interphase  is  provided  by  the 
difference  in  the  loading  and  unloading  rates  and  the 
incomplete  unloading  which  leaves  0.1-0.15  D/Pd  re¬ 
maining,  i.e.  characteristic  feature  (Hi).  This  behavior  is 
not  consistent  with  diffusion  control.  The  observed 
asymmetry  between  the  loading  and  unloading  time 
dependences  would  not  be  predicted  by  a  linear  diffu¬ 
sion  model.  Such  a  model  would  be  analogous  to 
ordinary  diffusion  with  constant  initial  bulk  concentra¬ 
tion,  yielding  a  solution  that  is  the  product  of  a  time- 
dependent  factor  containing  geometrical  information, 
but  not  initial  conditions,  and  a  function  of  the  differ¬ 
ence  between  the  initial  bulk  concentration  and  the 
surface  concentration  (see  e.g.  eqn.  2.4  1.  of  ref.  3,  p. 
97).  Although  this  asymmetry  would  be  consistent  with 
surface  control  as  predicted  by  the  model  if  e.g.  /c_4  » 
k^  in  eqn.  (IV)  of  ref.  1,  this  choice  does  not  reproduce 
well  the  other  characteristic  features  (i)  and  (ii)  and 
fails  to  predict  incomplete  unloading.  Further,  this 
asymmetry  cannot  be  explained  by  diffusion  alone  even 
if  <  Dp  by  about  an  order  of  magnitude  [4]  because 
of  the  required  continuity  of  flux  across  the  a  -*  /3 
dividing  plane.  The  continuity  of  flux  and  the  low 
D/Pd  ratio  for  the  a  phase  would  result  in  a  very  thin 
a  region  that  can  be  identified  with  the  interphase 
itself.  Consequently,  another  meehanism(s)  must  be 
considered  in  seeking  the  explanation  for  this  observa¬ 
tion.  For  example,  the  incomplete  unloading  at  a  CD 
of  60  mA  cm  ■  or  larger  may  be  due  to  blocking  of  the 
electrode  surface  by  an  oxide  layer  in  as  much  as  the 
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applied  current  would  exceed  the  diffusion  current 
associated  with  the  transport  of  absorbed  deuterium 
out  of  the  Pd  lattice. 

Evidently  there  exists  a  mechanismfs)  that  causes 
the  absorbed  deuterium  to  affect  the  loading  and  un¬ 
loading  rales  differently.  If  we  assign  an  active  role  to 
the  interphase  as  suggested  by  the  surface  morphology 
(see  the  subsequent  discussion  of  the  interphase),  then 
a  mechanism  is  provided.  To  illustrate  this,  the  electro¬ 
chemical  potential  of  the  adsorbed  deuterium  contains 
a  contribution  resulting  from  the  electrode  overpoten¬ 
tial  while  that  in  the  bulk  is  affected  by  contributions 
involving  the  electronic  as  well  as  mechanical  interac¬ 
tions.  Upon  current  reversal  (during  unloading)  the 
magnitude  of  the  driving  force  may  be  different  and 
the  asymmetry  would  fallow  naturally. 

2.3.  Effect  of  electrode  radius 

The  effect  of  the  electrode  radius  on  the  calculated 
charging  curves  is  shown  in  Fig.  4,  where  the  diffusion 
coefficient  for  0-PdD  was  used.  It  illustrates  agree¬ 
ment  with  the  experimental  data  (see  ref.  2,  Fig.  6) 
except  for  the  electrode  with  a  radius  of  0.1  cm.  The 
anomaly  in  the  latter  concerns  the  initial  portions  of 
the  measured  charging  curves  for  which  there  is  no 
observed  difference  in  slope  between  the  data  for  radii 
of  O.OS  and  0.1  cm,  requiring  that  the  initial  charging 
efficiency  be  proportional  to  the  electrode  radius  —  an 
unrealistic  condition.  Otherwise,  eqn.  (1),  which  pre¬ 
dicts  the  slopes  of  the  charging  curves  to  have  a  1/r 
dependeiKe  initially,  is  well  satisfied.  All  three  calcu¬ 
lated  curves  erdiibit  an  initial  charging  efficiency  of 
0.78,  whereas  the  measured  curves  of  ref.  2,  Fig.  6 
show  an  initial  charging  efficiency  of  about  0.8  for  the 
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Fig.  4.  Model  cilcalertoiit  of  loadinf  cones  as  a  function  of  ekc- 
trode  ladtaK  a,  OiK  cm;  b,  0.1  cm;  c,  0.2  cm;  current  density  120  mA 
cra~^;4nMion  coefficient  1.6x10'*  cm^  s~';  other  panuneten  ate 
listed  Hi  Table  1. 


Time/min 

Fig.  S.  Model  calculatioos  of  loading  curves  for  current  densities 
>/mA  cm"^  and  diffusion  coefficients  D/cm^  s"':  a,  ; « 240, 
D-  1.6x10'*;  b,y  - 240,  £)- 1.6x10'*:  c.y- 240,  0-1.6x10"’: 
d,;-15.  0-1.6x10"*;  e,y-15,  D  =  1.6X  10'*;  f,  y  -  15,  0-1.6 
X  10  ’;  electrode  radius  O.OS  cm;  other  parameters  are  listed  in 
Table  1. 

O.OS  and  0.2  cm  radii  data  and  an  unrealistically  high 
value  of  about  1.6  for  the  0.1  cm  radius  data.  Because 
the  reaction  rates  have  no  dependence  on  the  elec¬ 
trode  radius.  The  ntode'  predicts  no  dependence  of  the 
asymptotic  charging  level  on  the  radius  of  the  electrode 
and  apparently  none  is  observed. 

2.4.  Electrode-loading  mode 

Participation  of  the  interphase  in  the  course  of  the 
electrode  charging  can  be  examined  numerically  by 
changing  the  value  of  the  diffusion  coefficient  used  in 
the  model  computations.  As  the  diffiision  coefficient  is 
made  smaller,  a  shift  to  transport  control  must  occur. 
The  effect  of  varying  the  diffiision  coefficient  on  the 
charging  curves  is  illustrated  in  Fig.  5.  At  the  highest 
CD  a  significant  change  is  noted  in  the  shape  of  the 
charging  curve  for  the  diffusion  coefficient  of  1.6  x 
10"’  cm’  s"‘.  In  contrast,  at  the  lower  CDs  only  a 
small  change  is  achieved  for  the  same  changes  in  the 
diffusion  coefficient,  indicating  a  relative  increase  in 
surface  control  at  lower  CD  as  expected.  This  further 
indicates  that  a  transition  between  diffusion  and  sur¬ 
face  control  occurs  within  the  practical  range  of  cur¬ 
rent  densities  at  the  diffusion  coefficient  qperable  in 
the  system.  A  comparison  with  Fig.  5  of  ref.  2  suggests 
that  the  observed  charging  rates  were  limited  by  diffu- 
si(Mi  at  CDs  of  60  mA  cm"’  and  greater,  while  for 
lower  CDs  the  limitation  may  have  been  the  surface 
reactions. 

Figure  5  also  indicates  that  a  more  realistic  treat¬ 
ment  of  the  bulk  diffusion  to  include  differences  in 
diffusion  coefficients  for  a-  and  j8-PdD  would  not 
predict  significantly  different  charging  curves  except  at 
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higher  CDs,  where  there  is  diffusion  control  of  the 
charging  process. 

3.  Evolution  of  the  PdD/DjO  interphase 

The  surface  inhomogeneity  with  regard  to  the  ab¬ 
sorption  of  hydrogen  isotopes  and  the  composition  of 
the  interphase,  as  revealed  by  X-ray  diffraction  spec¬ 
troscopy,  provide  further  evidence  for  the  active  partic¬ 
ipation  of  the  interphase  in  the  course  of  charging-dis¬ 
charging  processes.  The  experimentally  observed  dif¬ 
ference  in  the  loading  and  unloading  rates  and  the 
presence  of  residual  deuterium  in  the  Pd  electrode 
suggest  that  other  processes,  beyond  those  contained  in 
the  model,  play  a  significant  role. 

3. 1.  Surface  inhomogeneity 

The  active  participations  of  the  Pd  surface  manifests 
itself  as  inhomogeneous  changes  in  the  surface  mor¬ 
phology  of  Pd  electrodes  during  their  prolonged  expo¬ 
sure  to  evolving  hydrogen  and  deuterium  as  reported 
by  Rolison  et  al.  [5].  Another  way  to  demonstrate  the 
inhomogeneity  of  a  Pd  surface,  particularly  with  regard 
to  absorption,  is  to  view  it  using  Nomarski  optics  (6)  so 
that  regions  of  preferred  absorption  can  be  differenti¬ 
ated.  An  example  of  the  experimental  arrangement  is 
illustrated  in  Fig.  6.  A  Pd  wire,  heat  treated  at  the 
recrystallization  temperature  to  ensure  large  grains, 
was  mounted  in  an  epoxy  resin,  polished  (diamond 
paste  on  silk  followed  by  a  light  polish  with  0.05  /am 
A1;0,  on  microcloth)  and  etched  (anodic  etch  in  .'>(Ky 
HCI).  A  hole  was  drilled  through  which  a  Pt  counter- 
electrode  was  inserted  and  the  cell  was  closed  with  a 
thin  Mylar  film.  The  so-construeted  cell  was  connected 
to  a  power  source  and  placed  in  a  metallograph  (Leeo 
Neophot  21)  equipped  with  Nomarski  optics.  Images  of 
the  electrode  surface  were  eontinuitusly  reei)rded  tm 
videotape.  The  observed  changes  within  a  single  grain, 
which  is  outlined,  are  illustrated  m  Fig.  7.  It  is  unlikely 
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Fig.  7.  Preferred  locations  of  D.  penetration  into  a  single  Pd  gram  as 
obtained  by  electronic  subtraction  of  Nomarski  images.  Aficr  calhod- 
ically  charging  the  Pd  electrode  al  .S  mA  for  .VI  min.  an  image  was 
recorded.  This  image  was  digitized  and  the  image  of  the  grain  al  rest 
potential  was  subtracted  out.  The  subtracted  image  is  shown  with  the 
individual  grain  outlined.  The  magnilicalion  of  the  microscope  is  \ 
400. 


that  the  changes  observed  in  the  grain  can  be  ai- 
tributed  to  the  surface  or  near-surfacc  cr.ieking  ih.n 
has  been  found  after  long  charging  times  at  high  ciu- 
renl  densities  [7).  Rather,  these  changes  arc  attributed 
to  localized  volume  c.xtcnsion  and  indicate  th.it  even 
within  a  single  grain  there  are  preferred  sites  ol  .ib- 
sorption  —  a  conclusion  also  reached  bv  Rolison  ei  al 

[5]. 

3.2.  Su[)en  hargcil  region 

The  present  model  employs  a  simple,  passive  intei- 
phase.  In  aclualily.  as  indicated  by  Bucur  .ind  Boia  (.sj. 
the  H  atoms  accumulate  in  the  inlerphase.  i  e  a  super 
charged  region  exists.  To  explore  further  the  I’d  D-O 
interphase,  vve  examined  the  time  evolutuni  oi  the  m 
situ  X-ray  spectra  covering  a  rather  narrow  In  i.inge 
trom  .55°  to  55°.  For  this  purptise  a  cell,  shown  in  lag 
S.  connected  to  a  potentiostat  (FAR  model  ’(>.’)  w.is 
placed  in  an  .X-ray  diffraction  system  (Rigaku  Rl  2oo 
FI).  The  cell  body  and  stem,  on  to  which  a  I’d  foil  K  u 
clad  on  one  side  to  ensure  containment  ol  tlie  deu 
lerium)  was  attached,  were  constructed  out  of  ,i  ehemi 
cally  inert  material  ( I  P.X).  The  cell  was  piovided  wnh 
a  thin  polyethylene  wiiuiow.  which  minimized  se.iitei 
ing  of  the  X-rays,  and  a  small  hole  in  the  io(’  ol  ihe  cell 
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Fig.  8.  Electrolytic  cell  for  in  situ  X-ray  spectroscopy:  A,  polyethy¬ 
lene  window;  B,  electrolyte  cavity;  C  cell  body  (TPX);  D,  Pt  counter¬ 
electrode;  E,  Pd  foil;  F,  Pd  reference  electrode. 


body,  which  permitted  both  the  escape  of  gas  and  the 
addition  of  electrolyte  as  needed.  The  cell  was  assem¬ 
bled  in  a  dry  box. 

Figure  9  shows  the  diffraction  spectra  for  the  Pd  foil 
at  rest  potential  (spectrum  A)  and  at  subsequent  times 
and  potentials  (spectra  B-E).  At  rest  potential  the 
polycrystalline  Pd  foil  exhibits  two  peaks:  one  at  40.377® 
and  another  at  46.924*  due  to  reflections  from  the  200 
and  111  planes  respectively.  After  2  h  of  charging  at 
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Fig.  9.  Progress  in  electrode  ioadiiig  by  X-ray  spectra:  A  metallic  Pd; 
B  and  C,  transition  from  meudlk  to  /8-PdD  phase;  D  and  E,  8-PdD 
phase  (see  text  for  overpotentials  and  charging  tiroes). 


1.5  V  negative  overpotential,  during  which  the  elec¬ 
trode  appearance  had  not  changed,  new  peaks  at 
37.758®  and  45.025®  emerged,  indicating  formation  of 
the  ^-PdD  phase  (spectra  B  and  C).  Comparison  of  the 
ratios  of  the  111  peaks  to  the  200  peaks  for  the  Pd  and 
/3-PdD  phase  leads  to  the  conclusion  that  deuterium 
preferentially  enters  the  Pd  lattice  through  the  111 
sites.  This  conclusion  agrees  well  with  that  based  on 
theoretical  calculations  [9].  With  prolonged  charging, 
the  silvery  color  of  metallic  palladium  changed  to  black. 
After  24  h  of  charging,  the  recorded  spectrum  shows 
only  two  peaks  attributed  to  the  ^-PdD  phase  (spec¬ 
trum  E).  The  1 1 1  and  200  reflection  peaks  have  shifted 
to  lower  values,  indicating  further  expansion  of  the 
lattice.  With  continued  charging,  these  peaks  have 
broadened  (spectra  D  and  E)  but  their  peak  positions 
have  stayed  relatively  constant.  We  interpret  both  this 
broadening  and  the  shift  to  lower  2B  angles  as  indicat¬ 
ing  the  presence  of  a  supercharged  layer. 

4.  Strengths  and  weaknesses  of  the  model 

Models  are  used  both  to  simulate  a  complex  reality 
by  simpler  arrangements  and  to  display  the  essential 
features  of  a  physical  system.  Often  even  a  simplified 
model  clarifies  the  interpretation  of  ambiguous  data 
and  avoids  lengthy  calculations  that  might  be  irrelevant 
to  the  problem  considered.  In  the  present  case  even 
our  simple  model  indicates  that  the  observed  charging 
behavior  of  Pd  electrodes  cannot  be  explained  solely  by 
the  diffusion  characteristics  of  deuterium  in  the  bulk 
but  must  include  consideration  of  the  surface  reactions 
as  well.  This  is  especially  true  at  lower  CDs. 

We  have  applied  a  model  where  the  structure  of  the 
interphase  as  well  as  the  operating  driving  forces  are 
determined  by  participating  processes  occurring  in  both 
contacting  phases.  However,  we  have  treated  the  inter¬ 
phase  as  a  passive  element,  i.e.  in  effect  assuming  that 
the  surface  is  homogeneous  with  respect  to  the  chemi¬ 
cal  potential.  Such  an  assumption  is  not  realistic  (see 
Fig.  7),  although  it  might  be  justified  on  the  basis  of 
the  existence  of  the  supercharged  region.  Gearly  an 
assumption  of  uniform  surface  coverage  of  the  Pd  rod 
cannot  be  made.  Even  so,  the  model  does  remarkably 
well  in  predicting  a  saturation  of  the  initial  charging 
rate  and  a  non-monotonic  dependence  of  the  asymp¬ 
totic  charging  level  on  CD.  The  time  required  to  obtain 
a  full  charge  is  also  predictable  with  reasonable  model 
parameters.  This  model  provides  a  framework  within 
which  questions  such  as  the  relative  influence  of  the 
Volmer  and  Heyrovsky-Horiuti  charge  transfer  paths 
as  a  function  of  CD  can  be  considered. 

Although  the  model  does  predict  some  small  asym¬ 
metry  in  loading  and  unloading,  it  does  not  appear  able 
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to  piedict  as  large  an  effect  as  was  observed  and  it 
does  not  predict  the  incomplete  electrode  unloading  at 
all.  Evidently,  a  more  realistic  model  should  incorpo¬ 
rate  an  active  interphase  that  will  provide  mechanisms 
for  asynunetric  loading  and  also  explain  why  the  charg¬ 
ing  rate  and  level  can  differ  from  electrode  to  elec¬ 
trode. 
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